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Abstract 

^ : i 

CN . We consider the problem of steering control for the systems of one spin k particle 

and two interacting homonuclear spin ^ particles in an electro-magnetic field. The 
describing models are bilinear systems whose state varies on the Lie group of special 
r**"- . unitary matrices of dimensions two and four, respectively. By performing decom- 

positions of Lie groups, taking into account the describing equations at hand, we 
\q ' derive control laws to steer the state of the system to any desired final configuration. 

Explicit formulas are given for the parameters involved in the control algorithms. 
Moreover, the proposed algorithms allow for arbitrary bounds on the magnitude of 
the controls and for some flexibility in the specification of the final time which must 
be greater than a given value but otherwise arbitrary. 
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1 Introduction 



In recent years, there has been a great amount of interest in the study of control of 
quantum mechanical systems (see e.g. [||], 0, |23|, |33[). This has been motivated 



by recent advances in the area of nuclear magnetic resonance and laser spectroscopy 
which have rendered possible the introduction of active control at the atomic level. A 
major motivation to study control of quantum mechanical system is given by quantum 
computation @. The information in a quantum computer is encoded in the state of a 
quantum system which has to be manipulated in order to initialize the computer, perform 
logic operations and measure the result of the computation. Existing techniques for 
the manipulation of the state of quantum systems only allow to perform very simple 
operations. On the other hand, the introduction of a control theoretic point of view 
promises to greatly increase the number of operations that can be implemented with a 
quantum system as well as their accuracy. 
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The state of a general multilevel quantum system is described by a vector \ip(t) >, in 
a finite dimensional Hilbert space. At every time t, we have 



\m>=x(t)\m>- a) 

The operator X(t) is the evolution operator solution of the Schrodinger equation 

ihX = HX, (2) 

with initial condition equal to the identity, and H is the Hamiltonian operator. In many 
experimental situations, such as nuclear magnetic resonance spectroscopy, the Hamilto- 
nian H has the form H := H + Y^=i HiUi(t) where Ut(t) are externally applied electro- 
magnetic fields which play the role of controls. This simplified model assume no interac- 
tion with the environment other than through the controls. In this situation, equation 
(|2]) can be written as 

m 

X(t) = AX(t) + BiXtyuiit), X(0) = I nxn . (3) 
i=i 

In (|J), A, .E?i,..., B m , are matrices in su(n), where n is the number of levels of the system 
under consideration. The solution of X(t), varies on the Lie group of special unitary 
matrices of dimension n, SU(n). In this paper, we are interested in the control of X in 
(H). This is a way to obtain the control of the state \ip > in (H) and also to control the 
operation in ([I]) to be performed on the state of the quantum system. 

Decompositions of Lie groups, and in particular of the Lie group SU(n), have been 
recently used to prescribe control laws and to study the controllability properties for 
multilevel quantum systems described by equations such as (|3|) (see e.g. |14]], [fL5H , p4 



25], f29|). To illustrate the basic idea, let us consider the simplest case of a two level 



quantum system described by 

X = AX + BXu. (4) 

The matrices A and B are in the Lie algebra su(2) and the initial condition is assumed 
to be the identity. The matrices A and B are linearly independent and therefore they 
generate the whole Lie algebra su(2), since su(2) has no two dimensional subalgebras. [] 
DefineQ 

k := ,I <A J A >. (6) 
V < B,B > K ' 



1 If A and B are not linearly independent, the solution of (^), with initial condition equal to the 
identity, can be written as 

X(t)=e At+B fo^ T)d \ (5) 

and, to obtain a desired final configuration, we can just select a control function u to make 
At + B J Q u{t)(1t equal to one of the logarithms of the desired target Xf, if Xf is in the Lie subgroup of 
SU(2), described by {X £ SU(2)\X = e Bs , s e R }. 

2 The inner product < •, • > is defined as < A, B >:= Tr(AB*) and it is equal to the Killing form on 
su{2) (see e.g. 0). 
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The quantity A; is a measure of the 'control authority' of the given system. The matrices 
A + kB and A — kB are orthogonal and therefore, if T\ is the unitary matrix which 
diagonalizes A + kB, we have 

Ti(A + kB)T* = —iXS z , (7) 

and 

Ti(A - kB)T* = -i(aS y + bS x ), (8) 
for some parameters A > 0, a, b not both zero. S x , S y and S z are the Pauli matrices 

*=-§(! i). < 9 > 
s »--=K°i o)- < io) 

The matrices —iS x , —iS y , and — zS^ form an orthogonal basis in su(2) and we have the 
basic commutation relations 

[ iS x , iSy] := "iS Z) [ iSy, iS z ] := iS x , [ iS z , iS x ] := iSy. (12) 
Define the unitary matrix 

T 2 := ( a g+ib J • (13) 

We have 

T 2 T X (A + kB)T*T£ = -i\S z , (14) 



t 2 t\(a - kB)T*r; = -iV^T¥s y . (15) 

In geometric terms, the transformation T\ is a rotation which ensures that the z axis is 
aligned along the direction specified by A + kB while T 2 is a rotation about the z axis 
to align the direction of rotation specified by —i(aS y + bS x ) with the y direction. The 
change of coordinates T 2 T\ shows that it is always possible to assume that A + kB and 
A — kB are proportional to the Pauli matrices S z and S y , respectively. 

Now assume we want to find a control input steering to a final configuration Xf (after 
the change of coordinates). We can first express Xf using the Euler parametrization of 
matrices in 577(2) (see e.g. |3"0fl ) 



with a, 7 € [0, 4m], (3 G [0, 7r]. Once the parameters in (|iq) are known, it is immediate 
to find a control function steering the state of the system to the value Xf. We have the 
following result 
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Theorem 1 The control piecewise constant and equal to k in the interval [0, —k in 
the interval (Z, % + v J +b - 2 ] and equal to k in the interval (\ + v J +b2 , 2 + v J +f)2 + f], 
steers the state of the system (j^) (in the new coordinates) to Xf in p6[j. 



Results similar in spirit to Theorem 1 have appeared in 0] [25]. In particular the 
paper |25[ contains a factorization of the Lie group SU(2) which allows control with 
arbitrarily bounded power. The authors deal with the case where A and B are orthogonal 
and then generalize to a large class of system on SU(2). 

Consideration of a bounded amplitude control is quite natural in applications like nu- 
clear magnetic resonance where the term A in (^) corresponds to a large constant magnetic 
field (so that the value of k in (|4]) is typically large) while the time varying transverse 
magnetic field u has small amplitude. In the next section, using a Lie group decomposi- 
tion generalizing Euler decomposition (|16D we will give a control algorithm for any system 
on SU(2) of the form (|j) which allows for arbitrary bounds on the control magnitude. 
Sections 3 through 5 are devoted to the more complicated case of two interacting spin | 
particles. Emphasis is given to the homonuclear case, which for us means that the two 
particles have the same gyromagnetic ratios. We derive, in Section 3 the model of two 
interacting spin | particles and prove some properties of the Lie algebra su{n) which are 
relevant for the controllability of this system. The controllability of the system of two 
homonuclear spin | particles is dealt with in detail in Section 4, where we give explicit 
expressions for the reachable sets at every time. This treatment gives information about 
what states can be reached at any given time. In Section 5, we present an algorithm 
to steer exactly the state of this system to any prescribed final condition with arbitrary 
bounds on the control amplitude. Conclusions are presented in Section 6. 



2 Control of two Level quantum Systems with Arbi- 
trarily Bounded Control 

Consider once again system (f|) that we rewrite here 

X = AX + BXu, (17) 

with A and B in su{2) and linearly independent, and X(0) = I 2 x2- Assume that the 
amplitude of the control u is bounded by \u\ < M. We will consider a piecewise constant 
control which attains only the values ±M (bang-bang). In order to specify the switching 
times for the bang-bang control, we will consider a decomposition of the Lie group SU(2) 
in terms of the matrices 

Zi := A + MB, (18) 
Z 2 := A - MB. (19) 

In the following we assume that Z\$ are linearly independent matrices in su{2) and 
we consider the following parameters: 
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• Ai 5 2 > are defined as the magnitudes of the purely imaginary eigenvalues of Z^ 2 . 
The one parameter subgroups associated to Z 12 are periodic with period and 
Ai,2 is a measure of the speed at which one moves on the one parameter subgroup 
corresponding to Z\ and Z 2) respectively. 

• The parameter ip is the cosine of the angle between Z-y and Z 2l namely we define 

V> == <Z V Z2> r- (20) 

< Zl, Zl >2< Z 2 , Z 2 >2 

The parameter ^ is not changed by a change of coordinates (rotations of the ref- 
erence frame) nor by a scaling of Ai and/or A 2 . If if) — then Z\ and Z 2 are 
orthogonal. If and only if ip — 1, Zi and Z 2 are proportional to each other which 
we have excluded so that in general we have < < 1. 

Let 7i be the unitary matrix that diagonalizes Z\ namely 

T X Z{T{ = -12XA, (21) 



and define 
With T 2 given by 

we have 



T X Z 2 T{ := -i(aS y + bS x + cS z ). (22) 



T 2 '■— [ n a+»b J , (23) 



T 2 T x Z{T{Tl = -i2X 1 S z , (24) 

T 2 T X Z 2 T{T 2 = -icS z - iVa 2 + b 2 S y . (25) 
Therefore there exists a unitary matrix FT := T 2 T\ such that 

W^W* = -i2AiS 2 , WZ 2 W* = -i&S 2 - Wa 2 + b 2 S y := -iD. (26) 



The parameters Ai, A 2 , ip and the matrix W will all play a role in determining the 
generalized Euler parameters in the factorization of elements of SU (2) in terms of Z x and 
Z 2 described in the following theorem. 

Theorem 2. Consider an arbitrary (target) matrix in SU{2), Xf and let a G [0,47r] 7 
(5 G [0, 7r] and 7 G [0, 47r] &e the Euler parameters of the matrix WXfW* , namely 

WXfW* = e -iS z a e -iS y (S e -iS z y_ ^7) 

Choose a positive integer m such that 

cos 2 {^-) > i, 2 . (28) 
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Then Xf has the following factorization: 



X f = e Zl ^ (e Zltl e Z2t2 e Zlt3 ) m e Zl ^ , (29) 



with 

and setting 
or 

tfA 2 t 2 = f, 
if cf) > 0, and 
ift<0. 




r^(- 2 (^)-« (30) 

:= tan _1 (— -?/>tan(A 2 t 2 )), (31) 



7T 

-sign(i/;)-, (32) 



(33 > 



2tt + . > 

h = h = (34) 



The proof of the above decomposition is presented in ||. The described decomposition 
of the Lie group SU (2) and calculation of the parameters can be easily extended to the 
Lie group SO (3). The relevant calculations for the latter case are presented in Appendix 
B. 

F. Lowenthal |0| first showed that the Lie group 577(2) is uniformly generated [[§] |32| 



by any two linearly independent matrices in Z\ Z 2 . This means that every element Xf of 
SU(2) can be written as the finite product of alternate elements of the one dimensional 
subgroups corresponding to Z\ and Z 2 , namely 

X f = e Zltl e Z2 * 2 e Zl * 3 e Z2 * 4 • • • e Zlt % (35) 

for some parameters ti,t 2 , ...,t s > and that, although s depends on Xf, it is uniformly 
bounded over SU(2). The contribution of Theorem 2 above is that we provide explicit 
formulas for the parameters tj. The maximum number of factors 's' (maximum over all 
of SU(2)) is called order of generation and it is the minimum number of factors needed 
to express all of the elements of SU(2) as in flSHp. The order of generation depends on 
the angle if) defined as in fl20D , between Z\ and Z 2 . In particular, it is minimal and equal 
to three if Z\ and Z 2 are orthogonal. F. Lowenthal in [|T^] has derived a formula which 
relates ip to the order of generation s in (|35|). In particular, s = 3, if ip = and s = f + 2 
if 

C0S {1) < |^| < cos(j^), (36) 

with / > 2. 
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It is interesting to compare the number of factors for the factorization described in 
Theorem 2 and the minimum required according to Lowenthal formula (p53). The number 
of factor required in (|29|) is 2m + 1 where m is the smallest positive integer satisfying 
fl28|). The worst case is when /3 = tt, which shows that the minimum number of factors 
needed to express all the matrices Xf in 577(2) in the form ( p9|) is 2m + 1 where m is the 
smallest integer satisfying 

|V| < cos(^). (37) 

A comparison of (RTj) with (^) shows (identifying / + 1 = 2m) that the number of 
factors needed in our algorithm is exactly the one given by F. Lowenthal, namely the 
smallest possible, in the case when / in ( |3T)| ) is odd and larger by just one if / is even. In 
fact, the derivation of Theorem 2 was in the spirit of F. Lowenthal proof of his uniform 



generation result in [Tj| [|19|| . F. Lowenthal uses stereographic projections, translates the 



problem to the induced subgroup of the Moebius group and uses some of his previous 



results [0. Since we are interested in the determination of the parameters involved 
in the factorization (p9|), we worked directly on the Lie group SU{2) and derived the 
factorization using explicit expressions of the matrices involved. 

Given the decomposition in Theorem 2, it is now immediate to find a control algorithm 
to steer the state of system fll7|). We have the following: 



Theorem 3 Consider the piecewise constant control equal to M for a time and then 
equal to M, —M, M for times £3, £2 end t\, respectively, m times, and then equal to M 
for an interval of time of length This control steers the state X of ( [Hp from the 
identity to Xf in (\2Q). 

Notice that Theorem 3 is a direct generalization of Theorem 2 that can be obtained 
as a special case if tp — 0, m is chosen equal to 1 and one sets 2Ai = A, 2\ 2 = V a 2 + b 2 . 
Notice also that , although the algorithm allows to reach a given state Xf at a given time, 
say Tf, one can consider the actual final time as arbitrary as long as it is greater than 
a given value. This will be very important in the algorithms of Section 5 where we use 
this flexibility to obtain a given state at a given time in a suitable rotating frame. Let us 
illustrate this point by assuming (w.l.g.) Z\ = —iS z , with S z in ([TT]) • Then we have 



-iSyTT 



If we let the system evolve as e Zlt , for time t, then steer to e lSy7T in time Ti, then let the 
system evolve as e Zlt for time i and then steer to e~' lSy7T in time T 2 , we obtain from ( |38l) 
the identity matrix in time 2i + 7\ + T 2 . If we follow this procedure after having driven 
the state of the system to Xf in time Tf, we steer to Xf in time Tf + T\ + T2 + 2t, with i 
arbitrary. More examples of the use of this procedure will be given in Section 5. 

We conclude this section by discussing how the number of factors (switches) needed 
in the described algorithm, which is, as we discussed above, essentially the minimum 
number needed, depends on the amplitude of the control. Notice, from formula (pg|), that 
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the number of factors 2m+l increases with the value of By substituting Z\ = A+MB, 
Z 2 = A — MB in ip in fl2Tf), after some elementary manipulations, we obtain 

<A,b> W ( 39 ) 

<B,B>) 



k 2 -M 2 
^(P + M 2 ) 2 - 4M 2 ( 



The parameter k is defined in (||). A study of the function in (|39|) shows that M = k 
achieves the minimum value of number of switches. The function is decreasing in the 
interval [0, A;) and increasing in [k,+oo). The number of switches and the number of 
factors in (|35|) tends to infinity as M goes to zero as well as M goes to +oo. Therefore, 
in order to minimize the number of switches, if the control u is bounded by \u\ < M and 
M > k, it is convenient to use u = ±k rather than u = ±M. In other terms, there is a 
natural value for the control, given by k in (|B|) which is the best to use even though we 
are allowed higher amplitude controls, as far as the number of switches is concerned. 

3 Model of two interacting spin ~ particles 

We now turn to the more complicate case of the system of two interacting spin | particles 
used in NMR experiments. Recent literature we will refer to on this topic can be found in 



141 , |[26|| . The papers [^TJ, [22 introduce system-theoretic aspects of NMR spectroscopy. 



The Hamiltonian of a system of two interacting spin | particles which interact with 
each-other, and are immersed in a driving electro- magnetic field, is given by || [10 



Hit) := Y, Wvk + l2hk)u k (t) + ( 7l J l2 + l2 hz)u z + JIuhz- (40) 

k=x,y 



The first term on the right hand side of ([40]) represents the interaction of the two particles 
with the x and y component of the external magnetic field, u x and u y , which are allowed 
to vary with time. The second term represents the interaction with the z component of 
the field, u z , which is kept constant in Nuclear Magnetic Resonance experiments. The 
last term represents the interaction between the two particles which is modeled with a 
scalar Ising term. The constants 71 and 72 are the gyromagnetic ratios of particle 1 and 
2, respectively. The constant J 7^ is the coupling constant between the two particles. 
For k — x,y, z, we have 

I lk := a k <g> 1, (41) 

and 

I 2k := 1 ® a k , (42) 

where a k , k = x,y, z, are the components of the spin operator in the x, y, z direction and 
1 is the identity operator. Also, we use the notation 

hkhj ■= Ok ® crj, j, k = x,y, z. (43) 
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Schrodinger equation for the evolution operator X for this system is given byf] 



X = -iH{t)X, (44) 
where H{t) is given in (|40|) . We consider the basis | + + >, | + — >, | — + >, | > (spin 1 



up, spin 2 up and spin 1 up, spin 2 down, and so on), in the underlying four dimensional 



Hilbert space. In this basis, the matrix representatives of the tensor products in ((41|), 
(f4"2"D, ( |4*3D are the Kronecker products of the 2x2 matrix representatives of the operators 
that appear as factors, where the matrix representatives of u x ,y,z are given by the Pauli 
matrices (Q) ( |T0D (|TT|). We can write system (f40|), ( f4*4|) in the form 



X = AX + B x Xu x {t) + B y Xuy{t) + B z u z X, (45) 

with 

A ■= -iJI lz I 2z , B x>ytZ = -i{jiS XjVjZ ® 1 + 72 1 ® S^z), (46) 

and 1 represents the 2x2 identity matrix. The state X of fl45|) with initial condition equal 
to the identity, varies on the Lie group SU(A). It is well known |TB] |27| that every state 
in 577(4) can be reached from the identity by varying the (arbitrarily bounded) control 
functions u x (t), u y (t), if and only if the Lie algebra generated by {A + B z u z , B x , B y }, that 
we will denote here by £, is equal to su(4). In this case the system is said to be controllable 
(see also |], ||, |yj for more explicit controllability criteria for quantum systems). More 
in general, the set of states that can be obtained with arbitrarily bounded controls for 
system fl45|) is the connected subgroup of SU (4) corresponding to the Lie algebra C Two 
cases can be considered: the heteronuclear (71 7^ 72) and the homonuclear (71 = 72). We 
will prove some general properties concerning the Lie algebra structure of these systems 
and in the next two sections we will study in some detail the controllability properties 
and give control algorithms for the system of two homonuclear spins. 



In order to study the structure of the Lie Algebra underlying system ( f4"o]) , we will need 
the following general result on the structure of the Lie Algebra su(n), for general n > 2. 
The proof is based on the Cartan decomposition of the Lie group SU{n) and it is 
presented in 0. A review of the Cartan decomposition for general semisimple Lie groups 
can be found in Appendix A. 

Lemma 4 The subalgebra so(n) and every other element f G su(n), f ^ so(n), generate 
su(n). 

Proof. See §. 

For an heteronuclear system (71 7^ 72), it is easily seen that, with repeated Lie brackets 
of B x and B y , it is possible to generate all the elements of the form —ilik, —ihk in ( PD 
(f42"D . Therefore, the Lie algebra generated by B x and B y is given by su{2) x su{2). It is 
known that su{2) x su(2) is isomorphic to so(4), therefore, we would like to use Lemma 
4 to conclude that, no matter what the matrix modeling the interaction A ^ so(4) is, the 

3 We are setting the Planck constant fi = 1. 
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Lie algebra generated by A + B z u z , B x and B y is equal to su(4). This requires a little 
care if no further information is provided. Although su(2) x su(2) is isomorphic to so(4), 
it might 'sit' in su(4) in a different manner from so(4). It turns out that this is not the 
case since su(2) x su{2) is in fact conjugate to so(4) via an element of U(4). Therefore 
we can conclude with the following result which is independent on how we have modeled 
the interaction in the matrix A. 

Theorem 5 For every system of two interacting heteronuclear spins the Lie algebra C is 
equal to su(n) and the system is controllable. 

For homonuclear spins the Lie algebra £ associated to the model depends on the type 
of interaction we consider. For a general interaction of the form A = —i(ali x l2 x + bli y l2 y + 
chzhz) the Lie algebra C is isomorphic to u(2), if a = b = c ^ and to u(3) if at least two 
of the coefficients a, b, c are different. Let us consider in detail the case of Ising interaction 
where a = b = 0, c = J^0. 

The matrices B Xty>z are proportional to — i(S X} y sX ®l+l®S Xi y )Z ) and B x and B y generate 
a subalgebra isomorphic to su(2) and so(3). For reasons that will be clear shortly, we call 
this subalgebra /C. Now write A := —iJIizhz as 

A = A 1 + l -D, (47) 

with 

iJ 

A\ := — — {2l\ z l2z — hxhx — hyhy), (48) 

and 

D:=-U hkhk- (49) 

k=x,y,z 

The matrix -j commutes with A\ and B XtVtZ (and therefore with the Lie algebra generated 
by them). Repeated Lie brackets of A\ and B x ^ y ^ z generate matrices of the form —i{Sk ® 
Sj + Sj Cg> Sfc), j 7^ k, j, k G {x, y, z} and —i(S r (g> S r — S m <S> S m ), r,m e {x, y, z} that span 
a vector space that we denote by V . Since A\ G V, D e C. We define Q := K © V, and 
we have [K, K] C JC, [JC, V] C and [P, P] C /C. Therefore 

£ = spcmA (50) 

and the subalgebra Q has a Cartan decomposition as described in Appendix A. 

More information can be obtained on L if we perform a change of coordinates diago- 
nalizing D in (|49|)^. This transforms the matrix D into D = —ijdiag(—3, 1, 1, 1) and all 



l X -> TXT* with 
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the matrices in Q into matrices of the form 



with i? e sw(3). In particular, for matrices in JC, R E so(3), while for matrices in "P, 
R £ S, the vector space of 3 x 3, zero trace, symmetric, purely imaginary matrices. 
Therefore, is isomorphic to su(3) (notice the fact that all of su(3) is generated can be 
obtained as an application of Lemma 4) and y plays the role of 2/3x3- In summary, we 
have for L (the symbol ~ indicates Lie algebra isomorphism). 

C = span{^-} © K, © V « w(3) = span{il 3x3 } © su(3) = span{il 3x3 } © so(3) © 5. (52) 

In the following two sections, we focus on controllability analysis and control algo- 
rithms for the system of two homonuclear spin | particles. A number of algorithms based 
on Lie group decompositions can be found in the literature for the heteronuclear case (see 
e.g. [14j|, [24}] , ||25|| , p6|). Starting with a decomposition in factors of the target matrix, 



some algorithms use the so called hard pulses (namely very high amplitude controls) to 
obtain approximately state transfer within the Lie group corresponding to the Lie algebra 
generated by the B matrices in a system like (f45|). This along with zero pulses, which 
make the state of the system vary on the one dimensional subgroup corresponding to the 
matrix A, can be shown to obtain all the possible targets, if the system is underlying 
a Cartan decomposition. While this approach has been shown to lead to time optimal 
controls |T4| , it has been criticized because of the practical feasibility and possible side 



effects of the hard pulses. The paper [2(| contains a detailed discussion and a comparison 



between hard pulses and soft pulses control as well as an algorithm for the control of 
the system of two heteronuclear spin with arbitrarily bounded control. We shall present 
in Section 5, an algorithm that, without any approximation, steers the system of two 
homonuclear spin | particles to any configuration. This algorithm will be an application 
of the decomposition of SU(2) described in Theorem 2. 



4 Controllability analysis 

Using the Cartan decomposition for the Lie group su(3) above described, it is possible 
to obtain more information about the controllability of the system of two homonuclear 
spins. In particular it is possible to obtain explicit expressions for the reachable sets with 
piecewise continuous (but not a priori bounded) control. To do so, we will use the main 
result of IJJ which is called there the 'Time Optimal Tori Theorem'. In the following, we 



will denote by 7Z(t) the set of states that can be reached from the identity in time t and 
by n{< T) the set U < t < T ^(t). We shall denote by SU{2) 2 the Lie group corresponding 
to the Lie algebra (isomorphic to su{2) and so(3)) generated by {— i(Sk © 1 + 1 © Sk), k = 
(x,y,z)}, namely the Lie group of matrices of the form L © L, with L £ SU(2). This is 
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a not faithful (2 to 1) representation of SU(2) and therefore isomorphic to £0(3). This 
can also be seen by a change of coordinates as pointed out in the previous section. 

Theorem 7. For the system of two homonuclear spins, let D be the matrix defined in 
For every T > 0, 

closTZ{T) = e^ DT {X f e SU(4)\X f = K 1 e aiAl+a2A2+aaA3 K 2 }, (53) 
with K12 G SU{2) 2 , 0^2,3 > 0, ot\ + a 2 + a 3 = T, and 

A 1 :=A- l -D = -l(2I lz I 2z - I lx I 2x - I ly I 2y ), (54) 
iJ 

A 2 := — — {2IiyI 2y — l\ x l 2x — Ii z I 2z ), (55) 
iJ 

A3 := — — {2I\ x I 2x — l\ z l 2z — l\ y l 2y ). (56) 

Proof. See §. 

The eigenvalues of the matrices A\, A 2 and A 3 are given by 

1 1 

Ai = A 2 = --Ji, A 3 = 0, A 4 = ~Ji, (57) 
o 6 

therefore the functions e Aji , j = 1,2,3 are periodic with period ¥jr. This shows that if 

T > in (p3]) clos1Z(T) = e^~G, where G denotes the Lie group (isomorphic to SU(3)) 
corresponding to the Lie algebra Q isomorphic to su(3) generated by A\, B x , B y and B z . 
With these notations, we have 

Theorem 8. IfT> then 

DT 

1Z(T) = e~G (58) 



Proof. See §. 



5 Control algorithms for two homonuclear spins 

In this section we give a constructive control algorithm for the system of two homonuclear 
interacting spin | particles, which, in the spirit of Theorem 3, allows for arbitrary bounds 
on the control. This algorithm is based on a Lie group decomposition of the Lie group 
SU(3), and on an application of the decomposition result of Theorem 2. . 

In the following we will set 7 = 71 = 72- Consider system ( fl"5f) with ji = 72 and define 

U := e~^ Dt X. (59) 
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Since D commutes with all the matrices in equation ( |45|) we have that U satisfies the 
equation 

U = A X U + B x Uu x + ByUu y + B z Uu z , (60) 

with A\ given in (|54]) (sse Theorem 7). The matrices Ax, B x , B y and B z generate a Lie 
algebra Q isomorphic to su(3) and in particular they are underlying a Cartan decomposi- 
tion of Q, in that A\ G V and B x , B y and B z generate /C, and Q = JC Q)V in the notations 
of Appendix A. It is also useful to perform a change of coordinates as described in Section 
3 diagonalizing the matrix D in ( PD as 

D = —i^diag(—3, 1, 1, 1). (61) 

This transforms A into the matrix A = —ijdiag(—l, —1, 1, 1), and B x , B y and B z into 
the matrices 762,3, — 762,4 and 763,4, respectively, where Sj^, j < k, denotes the matrix 
G so(4) with all zero entries except the (j,k)-th and the (k,j)-th that are equal to 1 



and —1 respectively. In these coordinates, the matrix defined in is given by A\ 



—ijdiag(0, — |, |, |). This shows that the Lie group on which the state of the system fl60"D 
varies is given by matrices of the form 

S=(J ° G ), (62) 
with G G SU(3). Therefore we can assume that the target matrix Xf has the form 

DT f 

X f = e—S f , (63) 



where D is now in the new coordinates, Tf > and Sf has the form in (|62|). In the system 
of differential equations in the new coordinates 

X = AX + B x u x X + B y u y X + B z u z X, (64) 

it will be convenient to scale the time t by a factor ^ and redefine the control functions 
as u x — > 7jw a; (jt), u y — > — 7jWj / (j£), m 2 — >■ 77M2, which gives simpler expressions for the 
matrices A, B XtVjZ . In particular, we have that the matrix B x ^ y ^ z become 

B x,y,z = (a c ° ) ( 65 ) 

\U 0(i, 2 ),(l,3),(2,3) / 

where 5^,^ are the standard basis matrices in so(3). The matrix D becomes D = 
—i^diag(—3, 1, 1, 1) and the matrix e Dt is periodic with period 4^ and therefore Tf in 
(|63|) is defined up to multiples of 47r. The matrix Ax takes the simple form A\ := A — y = 
idiag(0,2, —1, —1). We shall refer to this time scale and control variables in the rest of 
our treatment. 
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We now assume that Tf and Sf are given and we show how to steer the state of the 
system to the desired final configuration ( |6"3"D with arbitrary bound on the control. We 
first remark that the problem amounts to the following one: 

Problem 1 Steer the state of the system 

S = AiS + B x Su x + BySuy, (66) 

from the identity to e~ BzUz ^ Tf+ni7T ^Sf in time Tf + Anir for sufficiently large n, using 
piecewise constant controls of the form u x = ; u y = const, or u x = const, u y = 0. 

To see this, define 

S := e HD-B.*.yt X , (67) 



and using (|64T) , S satisfies 

S = A 1 S + e- B ^ t B y e B ^ t Su y {t) + e- B ^ 1 B x e B * Uzt Su x {t). (68) 



In deriving (|68|) we have used the fact that S 2 commutes with A and -D commutes which 
each one of the matrices in equation flBl"|). Now a direct calculation shows that if u x (t) = 
kcos(u z t) and u y = —ksin(u z t) then the last two elements in the above equation give 
kB x S, and if u x (t) = kcos{u z t) and u y (t) = ksin{u z t) the last two terms give kB y S. 
Therefore, if one keeps switching between these two types of controls the right hand side 
of equation (|68|) is alternatively of the form A%S+ kB x S and A\S + kB y S, where we have 
denoted by k any constant. Therefore the choice of the control reduces to the choices of 
the constants k. These calculations could have been carried out in the original coordinates 
and are often referred to in the physics literature as 'considering the system in a rotating 



frame' p6| . We emphasize that the calculations in the case considered here do not involve 
any approximation. 

Now, if we are able to steer the state of the system fl66|) or equivalently the state of 



with the above described controls to e B ^ u ^( T f+* nn ) Sf in time Tf + Ann, it follows 



from ( |67D that the state of the original system (|64"D is driven to es^" 1 " 4 ^ Sf = e s f S 



because of the periodicity of e Dt . 
Solution to Problem 1 

Notice that, in the new coordinates, the first rows and columns of the matrices Ai, 
B x and B y are zeros, so that we can consider the system as a 3 x 3 one with state varying 
on SU(3). 

We show next, how to obtain every matrix of the form e Bxt = e Sl - 2t (notations are for 
3x3 matrices). Set u y = and write A\ = —F + Ai, with F = —idiag(^, |, —1) and 
Ai = —idiag(^-, §,0). Setting 

S:=e Ft S (69) 
and noticing that F commutes with both A\ and B x , we obtain 

S = A 1 S + B x Su x . (70) 
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e Bxt 



Now, since the third rows and columns of system (|70| ) are all zeros the system is essentially 
2x2 and the matrices are essentially in su(2). Therefore, we can apply Theorems 2 and 
3 to steer the state of this system from the identity to the matrix^ 

/ cos(t) sin(t) 0\ 

-sin(t) cos(t) . (71) 
V \) 

Let T x = T x (t), the time needed to steer to e Bxt with bound on the control, according to 
Theorems 2 and 3. Now, notice that e~ Bx 2 A\e Bx ^ = —A\ and therefore, for every r > 0, 
we have 

Therefore, we can steer to the identity in time 

T Id (r):=2r + T x (^) + T x ( 3 -^) (73) 

with r > arbitrary We would like to steer the state of ( ffO|) to e Bxt in time 4wr, for 
some integer n. We can do that by steering in time T x (t) to e Bxt and then by steering 
(from the identity) to the identity in time Tj d {r) in (|73|), choosing r so that 

T x (t)+T Id (r)=Ann. (74) 

Notice that, from the decomposition of Theorem 2, the time T x (t) is uniformly bounded 
(as < t < 2tt) by a value T x that depends on the bound on the control. T x can be 
easily calculated by knowing the order of generation in terms of the two generating one- 
dimensional subgroups and knowing their periods. Therefore, we can always assume that 
n is chosen so that we can make (|74[) satisfied for every t, < t < 2n, for appropriate r, 
namely 4nvr > T x + T x {\) + T x (&). We will assume this to be the case in the followine. 

Now, since e Ft is periodic with period 4n, by steering the state S of ( |70|) to e Bxt in 
time Ann, using (|69|), we have steered the state of the system (|66| ) to e Bxt in time 4n7r. 

In a completely analogous way, one can obtain matrices of the form e Byt = e Sl ' 3t in 
time 4n7r (in fact the procedure is completely equivalent and the value of n is the same 
in the two cases). Moreover, noticing that 

one can obtain every matrix of the form e Bzt in time 3 x Ann. 



At this point, we recall a decomposition of the matrices Uf G SU(3) presented in |p0 |. 
Every matrix Uf G SU(3) can be written as 

U f = D(ai, a 2 , a 3 )U 12 (0i, cri)U 13 (9 2 , cr 2 )U 23 (0 3 , cr 3 ) , (76) 



5 An alternative here could have been to apply the algorithm of p4| appropriately modified to allow 
bounds on magnitude rather than power. This is possible since in this case the two matrices A\ and B x 
are orthogonal. 
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where the matrices Uki{0, cr), with k < I are the matrices whose submatrix intersection of 
the k— th and l—th rows and columns is given by 



so that, for example 
The matrix -D(cti, a2, a 3 ) is of the form 



sin(9)e la cos{6) 
( cos(9 x ) -ain(0i)e -fal 



sm(^i)e^ 1 cos(0i) | . (78) 
V 





/ e iai 








D(a 1 ,a 2 ,a 3 ) := 







e ia 2 







\ 








e ia 3 



(79) 

with ai + a 2 + «3 = 0. Once Uf is given, the parameters 9j,aj,atj, j = 1,2,3, can be 
easily calculated by an analytic procedure described in p0[ . Therefore, specifying Uf is 
equivalent to specifying the parameters 9j, atj, Oj. 

Now notice that 

U 12 (9 U a,) = e AlZ ^e B * 01 e Al ^. (80) 

Therefore the matrix Uu(9i, a\) can be obtained (assume w.l.g. <j\ > 0) by letting the 
system go with u x = u y = 0, for time and then drive the system to e^ 6 * 1 as explained 
above, and then let the system go again for time 2tt — setting m x = u y = again. The 
same thing can be done for the terms Ui3(9 2 , 02) and ^23(^3, cr 3 ), in ([76]), with e S:/612 and 
e Sz 3 replacing e BxBl , respectively. Therefore the last three factors in ([7j]) can be written 
as 

D(- ai , -a 2 , -a 3 )U f = e Altl e B * t2 e Alt3 e B y u e Alt5 e B * ts e Altr . (81) 

Our previous discussion shows how to calculate the parameters tj > 0, j = 1,...,7 and 
how to drive the identity to each one of the factors in (|8lD . Recalling periodicity of the 
matrix e Alt and the times to drive to the matrices e Bxt2 , e Byti , e Bztfi , calculated in the 
above discussion, the transfer to a matrix of the form (Ell) takes at most time 



Ti = 4 x 2tt + 2 x Ann + 3 x 4mr, (82) 

where the three terms on the right hand side refer to the matrices of the form e At , e Bxt 
and e Byt , and e Bzt , respectively. Now, the matrix D(ai, a 2 , a 3 ), with a\ + a 2 + 03 = 0, 
can be obtained as 

D(ai, a 2 , a 3 ) = e B ^e Altl e' B ^e Alt2 , (83) 

if 2t 2 — t\ — <x\ and 2t\ —t 2 = a 2 . Therefore the overall transfer of the state of ( pB| ) from 
the identity to Uf takes at most time 

f = f 1 + 2 x Ann + 2 x 2vr, (84) 
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with T% defined in ( |82|) . 

Now notice that, for every r > 0, we have 

e B v % eAl T e -B v z e B x % e A lT = /3x3 _ (g5) 

Therefore the Identity matrix can be obtained in time 3r + 4 x Ann, for arbitrary r > 0. 
Now choose n so that 

T f + Ann > f + A x Ann, (86) 

and express the matrix e - B ^ u ^( T f+ inn )Sj as Uf in (ff^). Then, we have showed how in time 
at most T we can steer to this final condition. Let T (< T) the time actually used to 
steer to e ~ Bz " z ( T ^ +4n7r )S/, and let us define 

i := Tf + Ann — T — 4 x 4n7r > 0, (87) 

from (|55|). Then by steering to the identity in time i + 4 x 4mr, we have driven the state 
of the system (|66|) to e~' Bz '" z ( T -'' +4 ' 17r )S'/, in time exactly Tf + 4n7r, as desired. □ 



6 Conclusions 

We have presented control algorithms and controllability analysis for the systems of one 
and two homonuclear spin ~ particles in a magnetic fields. In particular explicit expres- 
sions for the reachable sets have been obtained. 

The control algorithms are based on a Lie group decomposition of the Lie group 
SU(2) for which we have obtained explicit expressions of the parameters involved. All 
the algorithms described allow for arbitrary bounds on the controls and provide explicit 
expressions for the parameters involved in the control design. Moreover, the algorithms 
allow some flexibility in the final time that can be fixed a priori as long as it is greater 
than a given value. This is due to the right invariance of the systems considered and to the 
fact that the identity matrix can be obtained in sufficiently large but otherwise arbitrary 
time. Questions of optimization of the time transfer have not been considered here and 
will be object of future research. Also, the explicit generation of control protocols which 
use the decompositions of SU (2) described here as a basic building block for different and 
more complicated systems will be investigated. 
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Appendix A: Review of Cartan decomposition of Lie 
groups 

We review, in this Appendix some basic facts about decompositions of Lie groups based on 
symmetric spaces, also known as Cartan decomposition, that we will use in the following. 



We refer to the texts JTTJ [12 for further details and generalizations as well as for some 



of the terminology that we use here. Applications to control theory are considered in 
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Consider a semisimple Lie algebra Q and the corresponding connected Lie group G. 
Assume that Q admits a decomposition as a vector space 

Q = K, © V, (88) 

where K is subalgebra of Q, namely [/C, K] C /C. Moreover, assume that the following 
commutation relations hold among the elements of K, and V, 

[K,V)CV, (89) 

[V,V]CJC. (90) 

Denote by K the connected Lie group corresponding to JC and assume it is compact. 
Denote by P the image of V under the exponential map. Under the above assumptions 
and definitions, 

G = PK, (91) 

namely every element of G can be written as the product of an element of P and an 
element of K. 

Every element in V belongs to a Cartan subalgebra (namely a maximal Abelian sub- 
algebra) A [U whose dimension is called the rank of the Lie group G. Moreover, any 
two Cartan subalgebras are conjugate via elements in K, namely, if A' is another Cartan 
subalgebra, then there exists an element ifeK such that KAK^ 1 e A', for each A G A. 
As a consequence, every element P in P, can be written as P = KAK" 1 , where K is an 
element of K and A is an element of A, the connected Lie subgroup of G associated to 
A. Therefore, one can write 

G = KAK, (92) 

namely, every element in G can be written as the product of an element in K, an element 
in A and an element in K, in that order. 



The simplest example of Cartan decomposition (|9~2"D is Euler decomposition for SU{2). 
Consider the Pauli matrices in @, (|Tof ) , (|ll|), divided by i, —iS x , —iS y , —iS z . These 
matrices satisfy the commutation relations (|12|)- Therefore, one can take, for example, 
/C := span{S y }, and V : = span{S x , S z }, and A := S z , although every other combination 
would be possible. 

The decomposition fl92]) can be continued by decomposing K in the same fashion as 
G. Continuing this way, one end up with an expression of every element of the Lie group 
G in terms of the product of elements belonging to one dimensional subgroups. 
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Appendix B: A decomposition of 50(3) 



In the following we shall call Shk, h < k the matrix G so(3) which has zeros everywhere 
except in the h, k-th (k,h-th) entry which is equal to 1 (—1). Given a matrix Zi, there 
exists a matrix T\ G SO (3) such that 

T x Z x Tl := X 1 S 12 , (93) 

A x > 0. This can be easily seen by choosing T\ := [vi, v 2 , v 3 ] T , with v 3 such that vJZi = 
and with norm equal to one and V\ and v 2 such that \v\ 1 v 2 ,v 3 \ form an orthonormal 
basis in R 3 . We also set 

T X Z 2 T? := aS 12 + bS 13 + cS 23 . (94) 
Choose now T 2 := e Sl26 with 9 such that bcos{9) + csin(6) = 0. Then we have 

T 2 T x Z x TlTl := \\S\ 2 , (95) 

T 2 T x Z 2 TlT? := aS 12 + dS 23 , (96) 

for some parameter d ^ 0. Therefore, we can always assume that, in appropriate co- 
ordinates, the matrices Z\ and Z 2 have the form Z\ := A 15*12 and Z 2 := a5i2 + dS 23 , 
respectively. Moreover we can divide Z\ by Ai > (this has the only effect that, in the 
matrices if the form e Zl *, t has to be scaled by a factor Ai) and analogously we can divide 
Z 2 (in the new coordinates in (|9~6|)) by d and therefore the parameter t in the subgroup 
e Z2t has to be scaled by a factor d. Define p := §, we can assume, without loss of gener- 
ality, that the matrices Z x and Z 2 are given by Z\ := 5 12 and Z 2 := p5 i2 + 5 2 3, and we 
shall do so in the following. Notice also that the above manipulations do not modify the 
value of the parameter if} in (|20f) which is given, in terms of p, by 



_P_ 
a/T+7 



To express a matrix Xf G 50(3) as 

X f = e Zltl e Z2i2 e Zl ' 3 • • • e Zlt % (98) 
we first recall (see e.g. [|3"0[j ) that we can express every matrix X/ G 50(3) as 

with a, 7 G [0, 2ir] and /3 G [0,7r]. This is the classical Euler decomposition of a rotation. 
The Euler parameters a, (3 and 7 can be easily calculated by using an analytic procedure 
described in |[28|| . Now, since we have set Zi = Si 2 , the problem is to express every matrix 
of the form Xf = e S23/3 as in ( j98|) . To do that, we show how to express every matrix of 
the form e 523 ™ as 
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for sufficiently large m. An explicit calculation gives 



1+p c ps p—pc 



V 

-so 



1J 7)* 



(101) 



. -ep+p — s c+p 2 . 
\ >r j2 ^ r j2 / 



where rj := a/1 + p 2 , s := sm(r/t 2 ), c := cos{r]t2). Now, if we choose m so that 



2^ 2 - 1 < cos(-), 
m 



we can choose £2 so that 



c + p 2 cos(Vl + p 2 t 2 ) + P 2 



1 + p 2 



m 



(102) 



(103) 



Assume that this is the case and consider e Zit2 in ( |101|) with this choice. For brevity, 
let us call ctij its i,j-th element. By choosing t x so that cos (£1)013 + sm(£i)a23 = 0, we 
can make the 1, 3-th entry of the matrix e Zltl e z ' 2t2 equal to zero. The 2, 3— entry of this 
matrix is equal to sm(£i)o 13 + cos(£i)o 2 3 and it is equal to sin(-) up to the sign that can 
be changed by replacing £1 with £1 + ir. Let us now denote by a^- again the entries of the 
new matrix e Zltl e Z2t2 . By choosing £3 so that ansin^) + oi 2 cos(£3) = 0, we can make 
the 1, 2-th entry of the matrix e^^e^ 2 * 2 ^ 1 * 3 equal to zero, without affecting the 1, 3— th, 
2, 3— th and 3, 3-th entries (this procedure is reminiscent of the standard calculation of 
Euler angles in p§|). The resulting matrix is exactly e 523 ^. Therefore, by considering m 
factors, a constructive decomposition of 5*0(3) has been obtained. 
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